Tunneling in a bosonized Fermi liquid. Release of February 1, 2008 



L. S. Levitov'^'^ and A. V. Shytov'' 
(a) Massachusetts Institute of Technology, 12-112, 77 Massachusetts Ave., Cambridge, MA 02139 
(h) L. D. Landau Institute for Theoretical Physics, 2, Kosygin st., Moscow, 117334, Russia 

Abstract 

We apply the bosonization technique to the problem of tunneling in a Fermi 
liquid, and present a semiclassical theory of tunneling rates. To test the 
method, we derive and evaluate an expression for the tunneling current in the 
problem of two parallel Luttinger liquids. Next we study tunneling between 
parallel two dimensional Fermi liquids. In the absence of a magnetic field, 
the conductance I /V has a resonance peak at zero bias. Our estimate for the 
peak broadening agrees with the many-body theory results. In a magnetic 
field perpendicular to the plane, we find a tunneling gap linear in the field, 
in agreement with the experiment. In this regime the tunneling is spatially 
coherent over a distance of the order of the cyclotron radius. The coher- 
ence can be probed by a magnetic field parallel to the plane, which gives an 
Aharonov-Bohm phase to the tunneling amplitude. As a result, the current 

is an oscillatory function of the field. 
PACS numbers: 73.20.Dx, 73.40.Gk 
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I. INTRODUCTION 



High mobility electron gas in GaAs quantum wells forms a very clean Fermi liquid. By 
varying electron density one can change the effective strength of the Coulomb interaction, 
and study the effects of electron correlation. Of special interest is the behaviour of this 
system in a magnetic field, because it gives a unique opportunity to study the properties of 
a Fermi liquid in a magnetic field, with virtually no disorder effects. 

A new technique recently developed by Eisenstein is a tunneling experiment where two 
electron gases are confined in parallel wells separated by a thin barrier through which elec- 
trons can tunnel from one well to the other Unlike a conventional tunneling mea- 
surement, where electrons tunnel through a link between two reservoirs or through a weak 
place in a barrier, this system has a barrier that is very uniform. Because of this, and also 
due to a very small amount of an in-plane disorder, the momentum of a tunneling elec- 
tron is conserved. This gives two conservation laws for the two components of momentum, 
which, together with the usual energy conservation, highly restrict the phase space of pos- 
sible final states of the tunneling. As a result, the tunneling is of a resonance character, 
and the tunneling current displays a sharp peak. At close electron densities in the wells, 
the tunneling conductance measured as function of the density mismatch, has a peak of 
almost Lorentzian shape, and of the width determined by momentum relaxation time. This 
behaviour is accounted for by a free electron gas model Jll] . 

A very interesting change in the tunneling current is observed when a magnetic field 
is applied perpendicular to the plane of the barrier. In the field, the resonance peak is 
shifted away from zero bias by a certain amount, and a bias region is formed in which 
tunneling conductivity is almost totally suppressed. The physics of this "tunneling gap" can 
be understood in terms of the work one has to do to pull an electron out of one Fermi liquid, 
and to inject it in the liquid on the other side of the barrier The role of the magnetic 
field is that it squeezes electron states and makes them localized in the lateral dimension, 
which increases the energy of the Coulomb interaction of an electron with its neighbors. 
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The bias at which there is enough energy to transfer an electron is bigger than the Coulomb 
energy, and determines the gap and the position of the peak. 

The gap increases as a function of the magnetic field: linearly at weak fields, when there 



are many occupied Landau levels and more slowly at higher fields [|,H,|TO|. Let us 

mention that a similar tunneling gap, linear in the field, was seen by Ashoori in a tunneling 
capacitance experiment 0. 

So far, theoretical work has dealt mostly with the gap in the high field limit, where the 



problem has been treated using many different methods |14-[15|, all arriving at the same 
estimate for the gap: A ^ e^/e/, where e is the dielectric constant, and / is a characteristic 
interparticle separation in the system. (In the Quantum Hall regime, / is of the order 
of magnetic length Ib = (hc/eBY^"^ , and has a square root dependence on the electron 
density.) In the low field limit, there is a calculation by Aleiner, Baranger, and Glazman 
||16| , p!7[| , which treats the problem by using a hydrodynamical picture of an ideal conducting 
liquid in a magnetic field, which is supposed to be valid at scales much bigger than the 
cyclotron radius. In this work, the gap is found to scale as a square of the magnetic field, 
which is at odds with the experiment |]5|,|^. We think that the reason for disagreement is 
that the scale relevant for the problem is set by the screening length, which at a weak field is 
much less than the cyclotron radius. Therefore, an appropriate model should include Fermi 
liquid interaction effects, and thus the proper technique must be the Fermi liquid kinetic 
equation, rather than the long wavelength hydrodynamics. 

A conventional technique that deals with the problem of tunneling in a Fermi system 
is the many-body perturbation theory, which expresses the tunneling current in terms of 
one-particle spectral weights, and then calculates the weights from the Green's function 
perturbation series ||18[. However, because the interaction effects in this problem are in 
some sense always strong, and also because in a weak field limit the problem becomes 
semiclassical, a many-body perturbation theory calculation, although possible, is not very 
illuminating. Another option, less conventional in solid state theory, would be to use a 
time-dependent Hartry-Fock mean field equation in imaginary time, which is known in the 



theory of nuclear matter to be capable of dealing with collective tunneling phenomena such 
as nuclear fission |T^. 

In this paper we present an approach which better captures the semiclassical character 
of the problem at weak fields, and treats it as a tunneling in a "bosonized" Fermi liquid. By 
that we mean that the Fermi liquid is described in terms of a Fermi surface fluctuating in 
space and in time. Of course, such an approach is equivalent to the many-body perturbation 
theory carried out in the long wavelength limit. Our motivation for employing it, apart from 
interest in exploring a new technique, is in its simplicity, and in the possibility of establishing 
a clear relation with the phenomenological Fermi liquid theory. 

To clarify the last point let us recall that, when dealing with a one particle tunneling 
problem, it is customary to consider an equivalent classical problem and to evaluate the 
tunneling rate exponent as a classical action on a path in imaginary time. The Fermi liquid 
bosonization provides a natural extension of this approach to the problem of tunneling into 
a Fermi system. One can write an action which generates the Landau equation of the 
phenomenological Fermi liquid theory, and then look for least action paths in imaginary 
time. In such an approach one has Landau quasiparticles, represented by deviations of the 
Fermi surface from the spherical geometry, moving in imaginary time. To determine the 
tunneling rate exponent, one has to find a least action "bounce" solution of the Landau 
equation, which describes a particle injected into the Fermi liquid, and then later removed, 
and also the Fermi liquid response to the presence of the particle. 

Bosonization of a Fermi system is a well known method in the one dimensional problem. 
It has been pointed out by Luther that one can generalize the method to higher dimensions 
pO| , and since recent work by Haldane p3[ great interest has been generated in this new 



technique [^,^. We begin with a short review of the method; in dimension one, and in 
many dimensions (Sec. II). In Sec. Ill, we discuss how the tunneling problem can be treated 
in this formalism. 

In Sec. IV we turn to the specific problem of tunneling between two parallel Fermi 
liquids. There are some unusual features of this problem, from a technical viewpoint, that 



we discuss by using an example of two parallel one-dimensional Luttinger liquids. Due to 
the momentum conservation, the tunneling is spatially coherent along the barrier. In the 
instanton technique language, this is accounted for by summing over all points where the 
particle enters the liquid, and where the particle leaves the liquid. Such a summation gives 
an "entropical" contribution, which properly accounts for the conservation of momentum. 
By calculating the tunneling current between Luttinger liquids, we demonstrate that such a 
method agrees with the result obtained from the known exact Green's functions. Another 
feature, explicit in the one dimensional problem, is the role of momentum non- conserving 
interaction. If the Luttinger liquid Hamiltonian has only a term corresponding to scattering 
forward, there is a divergence in the tunneling current at zero bias, accompanied by vanishing 
of the current at any finite bias, which is expected from what has been said above. Such a 
singular behaviour can be regularized either by introducing scattering by a disorder, or by 
fluctuations of the barrier width. 

The entropical contribution, as well as the regularization by a momentum non-conserving 
scattering, are characteristic for the mult i- dimensional problem as well. In Sec. V we gen- 
eralize the observations of Sec. IV to the two dimensional problem, and also demonstrate a 
relation to the classical orbits. It turns out that, although we study dynamics in imaginary 
time, the spatial structure of the instanton is characterized by classical orbits found in real 
time, which are straight lines in the absence of a magnetic field, and Larmor circles in a 
magnetic field. 

Our new results, all contained in Sees. V and VI, can be summarized as follows. In the 
absence of the field, and at matched electron densities in the wells, there is a peak of the 
tunneling conductance I /V near zero bias. The peak width is determined by momentum 
non-conserving scattering, and the shape of the peak is Lorentzian. The physics here is 
similar to that of the peak of the equilibrium conductivity J/V/^o studied as a function of 
the density mismatch in the wells |]Tl] . 

At finite fields, which are much weaker than the Quantum Hall fields, the peak is shifted 
to a finite voltage and a tunneling gap is formed which scales linearly with the field. To 



compare this result to the hydrodynamical calculation |T^, we identify two terms in our 
instanton solution, accounting for a quasiparticle contribution to the gap, and for a collective 
Fermi liquid response effect. The latter is identical to the hydrodynamical theory result, 
however, the former is missing in it. The linear field dependence agrees with the recent 
experiment However, the dependence of the gap on the electron density, predicted to 

be a square root, has not been seen. 

Also, we find that the tunneling is spatially coherent over a large distance set by the 
cyclotron orbit radius, and propose to probe this coherence by a magnetic field parallel 
to the barrier. The advantage is that such a field does not affect the in-plane motion, 
however, there is a dramatic effect on the tunneling due to an Aharonov-Bohm phase in 
the tunneling amplitude. The Aharonov-Bohm effect destroys coherence of tunneling in the 
lateral dimension, and leads to oscillations of the current taken as a function of the parallel 
field. In the weak field limit these oscillations are given by a square of a Bessel function. 



II. BOSONIZATION FORMALISM 
A. One-dimensional bosonization. 

Bosonization of one dimensional Fermi systems is a powerful tool for the investigation 
of properties of low energy modes. It relies on the representation of the fermionic operators 
in terms of bosonic ones. In the low energy limit the current and the density of fermions 
become the only relevant variables, which makes the bosonization possible. Here we briefly 
review this method for the case of spinless fermions. 

For making a relation to the higher-dimensional problem more transparent, instead of 
the density p and the current j, it is convenient to introduce the densities of right and left 
movers, p± = \{p ± j) = which satisfy commutation relations: 

[p+(a;),P+(a;')] = [p-{x), p-{x)] =i5'{x-x') ; 

[p+(x),p_(x')] = 0. (2.1) 
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Then, one introduces the field operators 

0±(x) = ±E-e-"l^l/'-^^>±(p). (2.2) 

Here a is the ultraviolet cutoff set by the bandwidth. Conversely, the representation for 
fermionic operators has the form: 

^± = ^^e^^'=^^+^'^±(^) . (2.3) 
v27ra 



The Hamiltonian of interacting fermions is 

H = I dx 



— IVp ipj.—;:, ip 



dx dx J 

+ lT.^abPa{x)pb{x)] , (2.4) 

a,b 

where a,b = +, — , and W is a 2 x 2 matrix of Fourier components of the interaction U{x): 

= U = f/(0), and = = U{2kp). After the bosonization transformation, the 

Hamiltonian turns into 

H = ^JdxY^ {dM {Sab + F,,) {d^ct>,) , (2.5) 

a,b 

where Fab = ^iD^^ab- (Here = 1/2tivf is one half of the one-dimensional density of 

states.) One can interpret the entries of the matrix F as amplitudes of forward (-F++, F ) 

and backward (F+_,F_+) scattering. It will be seen later that the amplitudes Fab are ID 
analogs of the Landau function -F(ni, n2) of a Fermi liquid. 

For considering tunneling in a Fermi liquid we need an action corresponding to Eq. (|2.5| ). 
We derive it from the Hamiltonian by following the method of Ref. From Eq. (|2.1| ) we 
get the commutators of 0±(x), and use them to write the action: 

S=-^jdxdt {dt<P+d^<P+ - dt<P-d^<P-) - J Hdt (2.6) 

Note that our Lagrangian describes two chiral field theories interacting through backscat- 
tering. The Lagrangian of a chiral field theory [^] has the form: L^h = dx'P {dt + d^) (p- The 
solution of its equation of motion is 4>{x,t) = 0i(x — t) + 02(t), the second term 02 being a 
"ghost," which does not contribute to any observable (e.g., to p ~ dx4>)- One can say that 
the appearance of the 02 is the charge one pays for the locality of the Lagrangian. 



B. Multi-dimensional case 



It is well known that for a Fermi-system in any space dimension, the relevant low energy 
excitations are quasiparticles with momenta near the Fermi-surface. The Landau theory 
of a Fermi liquid enables one to describe all low energy physics in terms of a distribution 
of quasiparticles near the Fermi surface, i. e., in terms of the fluctuations of the Fermi 
surface geometry. Deformations of the Fermi surface are bosonic excitations, and one may 
try to bosonize Fermi-liquid in D > 1. It has been done in Refs. ||20| , |23| -p5|] by a procedure 
analogous to that used in dimension one. One introduces normal displacement Un(r) of 
the Fermi surface at the space point r, in the direction n, which means that "local Fermi 
momentum" is now n{pF + Un{r)). Displacements 'Un(r) represent a multidimensional analog 
of p±. Next, one defines 0n(r): 

w„(r) = (nV)0„(r) , (2.7) 
(cf. Eq.( p.2|) ). Then, a canonical operator of an electron can be written as 



iP = —L= V e^P^"'-+^'^(»''-) (2.8) 
v27ra n 

(cf. Eg. ( p. 3D ). Commutation rules for Mn(r), 

K(r),^/n'(r')] = ^5n,n'(nV)5(r-r') , (2.9) 

naturally lead to anticommuting ip^s. The Fermi-liquid Hamiltonian is expressed in terms 
of Un(r) as 



y"rfr'rfnun(r) ((1 + F)n)^(r) . (2.10) 
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Here dn is a normalized Fermi-surface area (in two dimensions, dn = d6/27T), and u is the 
density of states (z/2d = m/Trh^). The integral operator F describes quasiparticle interaction 
in the Fermi liquid. For simplicity, we assume the interaction to be of a density-density type: 
F = vU{v), where U is the interaction potential. It will be obvious to generalize the method 
to a more general interaction F{y — r', n — n'). 
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The action, by analogy with Eq.(|2.6|), is 

S = -^ J dtdrdn{nV)(l)nir)dt(l)nir) - J Hdt . (2.11) 

The equation of motion that follows from this action has the form of the Fermi liquid kinetic 
equation: 

[dt + VF{nV){l + F))unir) = . (2.12) 
It is not entirely trivial to include a magnetic field in this picture |23|. It cannot be done 



by a naive replacement V — V — |^A, because in the Fermi surface description one has 
to use the kinetic momentum mv, not the canonical one p = mv + eA/c. Hence, a gauge 
transformation of (p is not just (t>'y^{r) = 0n(r) + x(r)- To introduce a magnetic field, we 
consider the kinetic equation 

[dt + vpD {l + f)) u^{r) = , (2.13) 

where 

L) = n-V--B-nxVp. (2.14) 
c 



Following Ref. we call this expression a "covariant derivative," since it takes into account 
the curvature of the trajectory. If s is the length along the trajectory, then the rate of change 
of Mn(r) along it is given by 6u/6s = Du. 

Now we redefine so that u = Dcf) and rewrite the action to fit the correct equation of 
motion. 



S= '^"^ 



J drdndtD4>(dt + VF{l + F)D'^(j). (2.15) 



2 

This equation completes the description of the bosonized theory. 

III. TUNNELING FORMALISM 



Let us begin with some details on the system in which tunneling experiments are done. 
Typically, it consists of two quantum wells separated by a barrier of thickness d ~ 8 — 20nm. 
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The wells' width w is somewhat bigger: w ~ 15 — 30nm. The density of electrons in the wells 
is of the order of lO^^cm"^. Usually electrons fill only the lowest energy mode of transverse 
quantization. Because of this one can ignore the motion perpendicular to the plane, and 
treat electron dynamics as two-dimensional. The Coulomb screening length estimated 
from the density is of the order of 5 — lOnm. However, due to the finite well width w, 
effective interaction between electrons at distances shorter than w is reduced compared to 
the Coulomb potential, and this leads to the increase of the screening length up to a few w's. 
A magnetic field's strength varies from 0.5T to 5T and higher, corresponding to the filling 
factor varying from z/ ^ 10 to z/ ~ 1. For such fields, the magnetic length Ib = (hB/ecY^'^ 
is of the order of 10 — SOnm, which is less or comparable to the screening length. 

In this paper we are going to deal with cooperative effects that may strongly suppress 
charge tunneling. During the tunneling, the Fermi system has to accomodate the presence 
of a new electron, and the action of this process will completely control the tunneling rate. 
The single electron tunneling, which is a part of this more complex process, plays the role 
of an instant shake-up, followed by a slower relaxation in the many-body system. 

In such a problem, one can use the tunneling Hamiltonian formalism . By neglecting 
complications such as a finite well width, one can describe electrons in each well by a two- 
dimensional Fermi liquid Hamiltonian. The wells are coupled due to tunneling through the 
barrier, which is accounted for by a tunneling term Ht in the system Hamiltonian: 



Here to(r) is a tunneling matrix element. In a real system the barrier may be non-uniform, 
in which case to acquires a spatial dependence. 

The tunneling current due to a voltage drop across the barrier V can be found by a 
standard formalism ||18|. First, one may gauge the voltage away by transforming ipji — *■ 
ipRC^^^^. Then, one has to find the Matsubara transition amplitude. 




(3.1) 



K{v,t) = (T,7/^fl(r,2r)Vij(0,0)^L(r,^r)^i(0,0)). 



(3.2) 
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To extract real-time information from the Matsubara function one continues K to real 
time from the upper and lower half-planes. (Following Ref. we denote these functions 
as i^^*^^^(r, r)). The tunneling current / is determined by 

/ = 2 Re j dtdr dr' to(r)t*(r')$r-r',te-^"^* , (3.3) 

where $r,t = K>{Y,t) - K<{r,t). 

Now, let us derive an expression for the correlator in Eq.( |3.2| ). We use Eq.( p.8| ) to 
represent ip's. To evaluate the average, one has to write down a path integral of the product 
of four exponentials with the weight exp(— 5). (From now on the time is imaginary.) This 
Gaussian integral can be evaluated by a saddle-point method, so that 

K = -i- Y: e-^^""'"^^ , (3.4) 
(2vra)2 

where S{no, ni) is the saddle-point action 

5(no, ni) = '^{JD-' {tdt + vp{l + F)d)"' J) . (3.5) 
(...) denotes / . . . dndrdt, and 

Jr^l) = ± ((5(r)5(t)5„,„„ - S{r - R)5(t - r)5„,„J (3.6) 

Eq.( |3.5|) has a simple interpretation in terms of the kinetic equation ( p.l3| ). Suppose an 
electron is instantly transferred across the barrier at r = 0, t = 0, and then is transferred 
back at r = R, t = r. To incorporate the electron transfer, one must add a source term to 
the right-hand side of the kinetic equation ( p.l3| ): 

^Udt + D{l + F))u = iJ. (3.7) 

Then, one notes that Eq.( |3.5| ) simply gives the action for this solution, describing the 
Fermi liquid accomodation as an injection and removal of an electron. Finally, accord- 
ing to Eq. ( p.3| ), one has to integrate this action over R and r, with the weight e~*^^*. This 
weight can be interpreted as an additional term in the action, which accounts for the change 
in the potential energy during the tunneling time t due to the voltage source. 
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In this way, one may interpret Eq.( |3.5| ) as the action for a bouncing trajectory ||2^, and 
treat Eq.(p.3|) as a sum of the amphtudes of all the different bouncing trajectories. Quite 
often this summation can be done simply by taking the least action trajectory. However, 
in the case of tunneling between parallel layers this conventional procedure would lead to 
an incorrect result, because there exists a great number of trajectories and their "entropy" 
contribution becomes essential. As a technical remark, we will see that the summation over 
initial and final points is responsible for momentum conservation and can drastically change 
the answer. 

In the next two sections, we apply Eq. (|3.5|) to the problem of tunneling in a one- 
dimensional Luttinger liquid and in a two-dimensional Fermi liquid. 



IV. ONE DIMENSIONAL PROBLEM 



A. Clean system 



To illustrate how the formalism works, we first consider tunneling in a one-dimensional 
Fermi-system. Let us take two parallel Luttinger liquids with a tunneling coupling be- 
tween them. For simplicity we include only interaction inside the wire, but generalization 
is straightforward. Then Eqs.(X3) and ( ^ ) simplify to (we introduce matrix notation, 
corresponding to left and right movers) K 



^2 and 



S = 7t{JD-^ (idt + vf{1 + F)d) J) 



(4.1) 



where 



Jr(L) = ± imHx) - S(t - t)S{x - X)) 



(A 



vV 



(4.2) 





il 1 




Fi F2 






Fr{l) = 




^0-1 J 




^ F2 Fi ^ 



(4.3) 
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To find 5* in Eq.( |4.1| ), one introduces Fourier transform and carries out matrix inversion. It 
gives 

, r dkdu to — ikvpil + Fi) . r, cut + kX . ^, 

Here v = vf\J (1 + -^i)^ — F2 is tlie velocity of excitations in tlie Luttinger liquid ||2^. After 
integration over uj and fc, and summation over left and right movers (and cutting off the 
divergent integral at 1/a), one arrives at 

^2{/3-2) ^2 _ ^2^2 
K{x,t) = / 2 ^ 2 2^/3 ' (^•^) 



One can compare this expression to the product of two exact Green's functions of the 
Luttinger liquid |^|. Note that for a chiral liquid (i.e., with no backscattering, F2 = 0) one 



has /3 = 2, and only Fermi-velocity is renormalized. Note also that (3 does not depend on 
the sign of Fi and F2 to the first nonvanishing order of the perturbation theory in Fi and 
F2. This happens because of an approximate "particle-hole symmetry": if the interaction is 
attractive, the work should be done to remove a particle, while for a repulsive interaction 
the same work is required to inject a particle. 

To get the tunneling current, one has to compute the integral in Eg. ( p. 31) . On dimensional 
grounds, one expects the current to behave as a power law: / ~ \/2(/3~2) Explicit calculation 
(see Appendix A) gives 

m = C{P) ^ (^) , (4.7) 

1 (/3 - 2) r (/3 - f ) 

where T{f3) is the F-function. However, in an apparent contradiction to the dimensional 



estimate, for the simplest case of free fermions, Eq. (^l7| ) gives (3 = 2, C{(3) = 0, and I{V) = 0. 



The same cancellation of current also occurs in the absence of backscattering. 
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This cancellation has a clear physical reason. It has been pointed out p[ that, in a gas 
of free fermions the energy and momentum conservation prohibit tunneling at any nonzero 
voltage. This argument works even for an interacting system in the absence of backscatter- 
ing. However, it fails when the backscattering is present, because of the possibility of an 
electron scattering to another branch of the spectrum, with the 2pF change in momentum. 

Therefore, to get a finite tunneling current for a system with no backscattering, one 
needs to introduce some mechanism which changes the momentum of the tunneling electron 
(i.e., breaks translation invariance). There are two mechanisms in real systems: 

(i) scattering by a disorder; 

(ii) spatially fluctuating barrier width; 

or, more formally, a spatial dependence of the tunneling matrix element. If the barrier is not 
uniform, the tunneling occurs preferentially where it is low or narrow. Below we consider 
both mechanisms. 

B. Scattering by a disorder 

We consider here only weak disorder, so that electron motion can be treated as ballistic. 
This is the case when the characteristic time of motion is much smaller than the scattering 
time. In this limit one may take scattering into account semiclassically. Namely, since the 
probability for the electron to move by a distance s without being scattered is equal to 
exp{—s/ls) {Is is scattering length), we can simply multiply the amplitude K by this factor. 
After such a "regularization by scattering," one gets 

^rix, t) = $(a;, t) exp(-|x|//,). (4.9) 

One can compare this result to the perturbation theory expression for the Green's function 
in a random potential. In a ladder approximation the Green's function average over 
disorder differs from that in a clean system by the factor exp(— r/2Zs). 

After substituting Eq.(f4.9|) in Eq.(|3.3|), and integrating over x and t, one arrives at 
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2e|toP eV 

7ch {evy + {n/Ts) 



2 ' 



(4.10) 



where = Ig/vp- So, the effect of elastic scattering is that the conductance I/V becomes a 
Lorentzian function. 

Note that at small V Eq. f^loD gives Ohm's law, which differs from the power law ( [1.71) . 
This occurs because of the presence of an additional dimensional parameter Ig. 

Note also that the tunneling conductance in the Ohmic regime is proportional to Ig- Such 
a result is unusual for a tunneling problem. For example, for tunneling in a point contact, the 
I-V dependence is also Ohmic, but the tunneling conductance is expressed solely through the 
density of states, and the mean free scattering path Ig does not enter. Here, the appearance 
of Ig reflects an important feature of tunneling between two layers: together, the energy and 
momentum conservation prohibit tunneling at any nonzero voltage, and due only to elastic 
scattering does the tunneling become possible. One also may note that the behaviour of the 
tunneling current as ~* oo is extremely singular. 

To summarize, the disorder leads to the broadening of the resonance tunneling peak at 
= into a Lorentzian peak of the width h/Tg. 



Since the probability to tunnel depends on the width of the barrier exponentially, the 
tunneling through a barrier of a fluctuating thickness occurs at the points where the barrier 
is the thinnest. Near such points one can use the following model for the tunneling matrix 
element: 



Here a is the lateral scale of the barrier fluctuation. Since the motion of an electron in this 
model is always ballistic, one may assume an arbitrary voltage. By doing the integral in 
Eq.( p.3|) , one gets 



C. A non-uniform barrier 




(4.11) 
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/W = ~eVexp('-f^)V (4,12) 



TT V i 



Note that, compared to the previous case, the suppression of tunnehng at large V is stronger. 
Also, in the limit a — > 0, the barrier turns into a point contact, and one recovers from 
Eq.( 4.12 ) the usual Ohm's law. 

To summarize, the semiclassical treatment of tunneling through an extended barrier is 
different from that in a point contact. We will encounter a similar situation in the discussion 
of a two dimensional problem, so, to facilitate a comparison, let us list the main things we 
learn from the one dimensional tunneling problem: 

(i) The summation over initial and final points is an essential feature of the formalism, and 
can change the answer drastically; 

(ii) The tunneling current for free particles has to be regularized by breaking translational 
invariance; and 

(iii) The answer is sensitive to the regularization scheme, which has to be chosen on physical 
grounds. 



V. TWO DIMENSIONAL PROBLEM: NO INTERACTION 

It turns out that for 2D electrons in a magnetic field the technically difficult part is to 
carry out the semiclassical calculation for free particles. After that, taking into account the 
interaction will be straightforward. 



A. Relation of the least action paths to classical orbits 

Now, we will consider free electrons, F = 0, and evaluate the action ( p.5|) . It is worthwhile 
to mention that because of the semiclassical character of the problem, the calculation will 
naturally lead to a relation with the cyclotron orbits of the classical problem. To deal with 
the operator in ( p.5| ), it is useful to use Fourier-transform 
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= E'^fc-e-^"*^*i''' (5.1) 

k,ijj 

and to look for a solution to the equation 

itdt + D)G{e,e') = 5{9-e'), (5.2) 

where 6 is polar angle: n = cos 6*1 + sin^j. The solution is given by 

in{e-e')+ikx{n-n')Rc 

G(e,e') = y ^ — (5.3) 

where lOc = eB/mc is Larmor frequency and Rc = vp/ujc is Larmor radius. 
Now, it is straightforward to evaluate the saddle-point action: 

4 /■ (fkduo ^ 1 

Ofree = 777ITT ' ■ ^ 



v J (27r)3 inuJc{uj + inuJc 
sin^ ^(?2(6'o-6'i) + k(zx (no-ni)i?c+R) - icur) . (5.4) 



At this point we encounter a difficulty: the expression (|5.4|) is formally divergent. For 
electrons without a magnetic field it has been noted |^0|] that the action is finite only when 



there exists a classical trajectory with given momentum going from to R. If this is not 
the case, the transition amplitude becomes zero. Now, we will see that the same situation 
occurs in a magnetic field. 

To begin with, one divergence comes from the 0— th term of the sum in Eq. ( |5.4]) . This 
divergence will be eliminated, if we impose the condition 

z X (no - ni)i?c + R- = 0. (5.5) 



In the classical language, Eq.( ^.5|) means that the initial and final points belong to the same 
cyclotron orbit. However, this does not solve the whole problem, because now we have a 
diverging k-integration, which gives 

(5(z X (no - ni)i?c + i?) = 5(0) . (5.6) 

To regularize this (5-function, let us note that Eq. ( ^.5] ) has semiclassical accuracy, because it 
implies a definite trajectory of the particle. However, in a quantum problem there is a finite 
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thickness to anything, and thus the trajectory is washed out a bit, so that the (5-function 
effectively has a finite width. To determine this width, let us recall that, in a semiclassical 
picture, different states of a particle can be associated with the cells in the phase space of 
the volume {27ThY, where d is the space dimension. From this, a regularization of the phase 
space 5-function follows: 

Sp{0)5M = ^ , (5.7) 

where the subscripts denote the arguments of the 5-functions. To regularize the 5-function 
in Eg. ( ^.61) , we rewrite it in terms of the coordinate and momentum: 

Here, we drop higher order terms in the 5— function argument, since we are interested in 
dimension only. 

Therefore, we find that: 

(i) The action is finite only if there exists a corresponding classical trajectory; 

(ii) The regularization of the action can be done semiclassically. 

Let us mention a similarity of our regularization to the regularization of S^^\p) in the 
scattering amplitude in quantum field theory ^0|, done by a finite volume. (A finite volume 



problem with periodic boundary conditions, leads to a lattice in the momentum space, and 
a lattice 5-function is finite.) 

Now, the regularized action is finite: 

_ f duj 1 - cos {n{9o - Oi) - ujt) 

The integral over uo can be done by the residue method, and the sum can be calculated by 
using the identity 

^ = log^^, (Rez>0), (5.10) 

^1 n e^-l 

and cutting it off for z = Q divergent part at rimax = ^c/vpot (Recall that a is a short time 
cutoff of the order of inverse Fermi-energy or bandwidth). Finally, the answer reads 



K(R,t)= y ( ^ —-] , (5.11) 

where the sum is restricted to the points no, rii, R of the cyclotron orbit given by Eq. (|5.5| ). 

From Eg. ( [5 .31 ), a simple intuitive picture of tunneling follows. If one looks for a solution 
to the kinetic equation ( |3.7| ), one finds that it corresponds to the quasiparticle uniformly 
distributed along a cyclotron orbit at uJcT ^ 1. Clearly, it simply means that in a transla- 
tionally invariant system only the tunneling between cyclotron orbits with the same centers 
is allowed, and thus one can think of a tunneling electron spread over a cyclotron orbit and 
hopping to an identical orbit in the other plane. 



B. Semiclassical Jacobian 



The sum in Eq. ( ^.llD still requires some care, because of the semiclassical accuracy of 



Eq.([5.5|). Naively, one would simply take as a final answer the expression (|5.11|) for no and 



ni satisfying Eq.( |5.5| ), but this turns out to be incorrect. It is known |^ that a semiclassical 



propagator contains as a prefactor the Jacobian dpi/dqj, describing the divergence rate of 
two trajectories going from the same point with different momenta. The role of this Jacobian 
is analogous to the role of the prefactor of the semiclassical wave function in dimension one: 
tp exp(zS'). 

Below, we present an "intuitive" way to get the Jacobian from the summation over no, 
ni in Eq.(pll). 



To carry out the summation, let us consider the initial wave packet of the width Sx 
in a direction transverse to the trajectory. This packet has, therefore, a transversal mo- 
mentum distribution of the width 6p = 2TTh/6x. For free motion in a magnetic field, the 
corresponding uncertainty of the final point of the trajectory in the transverse direction is 
6x' = 6p/{pFRcSi'^{6i — ^o))- To have non-zero probability of finding the particle at R, 
6x' < 6x is required. By putting all together, one derives the semiclassical interval of angles 
over which the sum should be taken: 
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59 



5p 



(5.12) 



Pf 



PFRcSm{9i - Oq) 



From this, one can get the Green's function. By adding the amphtudes for (no,ni) and 
(—111, —no), one has 



where ^(r) = 6q — 6i = 2 sin ^r/2Rc. Now, for verification, one can check Eq.(|l|) with 
the semiclassical expansion of the exact Green's function. 

Finally, let us remark that the inability of our semiclassical method to automatically 
generate a correct Jacobian is a general property of a theory with linearized spectrum. The 
reason is that, in such a theory wavepackets do not spread, due to a linear dispersion relation, 
and thus the Jacobian must equal one. For example, in a one dimensional problem with 
a linear spectrum, in the semiclassical expression for the wavefunction the prefactor 
is absent. Taken literally, this leads to various anomalies, e.g., to the absence of a linear 
response to an electric field. In the one dimensional problem, this drawback is eliminated by 
incorporating a chiral anomaly into the theory. The bosonized theory reflects the anomaly 
as the Schwinger term in the right hand side of Eq. (|2.1| ). The multidimensional problem, 
however, does not incorporate it in a proper way. More exactly, only the "longitudinal" part 
of it is present in Eq. (|2.9| ), because the Fermi surface is assumed to be locally flat, while its 
actual non-zero curvature is responsible for the "transverse" part that leads to the Jacobian 
obtained above. Also, in our opinion, the same physical reason applies to all the divergences 
discussed in the previous subsection. 




(5.13) 



C. Regular izat ion by scattering 



We are going to introduce elastic scattering to regularize the expression for the tunneling 
current, similar to how it was done in one dimension. In a higher dimension, it turns out 
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to be a more complicated procedure, because of the uncertainty of the particle path going 
to a given final point r. In one dimension, the path length would be simply |r|, however, 
in two dimensions, since the motion in the magnetic field is periodic, the particle can make 
0, 1, 2, . . . complete circles before arriving at r. One can resolve this difficulty, by using the 
identity 

1 °° 1 

and expanding the amplitude K in the following series form: 



'S9 UJr. \ 



2 



T^TT^ X (5.15) 

By comparing it to the exact one-dimensional Green's function, one can see that the l—th 
term in this expression can be interpreted as an amplitude of / rotations before arriving at r. 
In this way, the path for each term is si — Rc\9 + 27rZ|. Therefore, due to elastic scattering, 
each term is reduced by exp(— s;//^). In this way, the regularized amplitude takes the form: 

Kr=(^^]\ (5.16) 

The scattering acts as a cutoff to this sum. Effectively, it eliminates the contributions of 
the trajectories which wind around the orbit so many times that their length exceeds the 
scattering length Ig. 

In order to be able to apply this method, one has to have a ballistic situation, Is- If 

this is not the case, the use of a coUisionless Fermi-liquid equation to determine the motion is 
incorrect. Therefore, this method enables one to deal only with a classically strong magnetic 
field. 
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D. Calculation of tunneling current 



Finite width of the experimentally observed peak has been attributed to the fluctuations 
of barrier width |^ . Here we calculate the effect on tunneling current due an in- plane elastic 
scattering. We expect that the broadening of the peak due to momentum non- conserving 
scattering (no matter what origin) will be similar to that found in Sec. IV for the one- 
dimensional problem. 

To get the tunneling current from Eqs.(^) and ( |5.16 ), we continue K to real time: 



K>(r,t) - i^<(r,t) = ^\ (^)' 

Pi^itc sm ty(rj \2tcvf/ 

Y.S\9~ u,t + 27r/)e-l^+2-'l^-/'^ (5. 17) 



Now one plugs this expression into Eq. (|3.3| ) and integrates it. Here it is useful to reduce the 
integration over r to the integration over 6: (fr = 2'Krdr = 27ri?, sin 6d6, since 2RcSm6/2 = 
r. The flnal answer reads 

where = h/vp- Note that this expression is identical, up to a prefactor, to the result 
for free electrons in the dimension one, given by Eq.( |4.1CI| ), and the same discussion applies 
here. 

VI. INTERACTING ELECTRONS IN A MAGNETIC FIELD 

A. Effective interaction 

In this subsection we treat F in Eq.( p.5|) . We start with an expansion of Eq.( p.5|) in 
powers of F, and rewrite it as a sum of two terms: 

1 11^1 , , 

+ -7-F^^-71^^ (6-1) 



D{LO + D{l+Fj) D{u + D) u + D uj+D' 

where 
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T = F + F^—F + --- = F + F^—T (6.2) 

The first term on the right hand side of Eq. ( |6.1| ) yields the free particle action ( pl9| ) studied 
above. The "effective" interaction T in the second term is defined by the Fermi liquid ladder, 
given by Eq.( |6.2| ). As usual, Eq.( |6.2| ) can be treated as a Dyson equation. It describes the 
screening of the interaction F by the electron gas. 
To solve Eq.( |6.2| ), we rewrite it as 

T{\\F) = F \F^^f (6.3) 
LiJ -\- L) 

and use Eq.(|5.3|) for (cu + D)^^ If F does not depend on angles (like for a density-density 
type interaction), a convolution gives 

(1 + F^)Tk^ = f 1 + 5: '^y.^^^ Tu}\ (6.4) 



n 



Here J„ are Bessel functions. 

Let us solve Eq.( |6.4| ) for small uj <^ uJc- The range of lo is set by relevant voltage, which 
in its turn is determined by the effective interaction. We will see below that the latter is 
small compared to TiiOci ^iid this implies that the voltage V is small: eV <^ huc- 

One has to take the limit uj carefully, because of the singular nature of the operator 
D^^. If one simply sets u = in Eq. (|6.2|) , one incorrectly gets = F/{l + F). However, the 
presence of a finite u changes the result. One may note that the limits n — and u in 
the 0— th term of the sum in Eq. (|6.4|) do not commute. The origin of such behaviour lies in 
the existence of an infinite number of stationary solutions to the kinetic equation Eg. ( p.l3| ). 
Indeed, if one neglects by the time derivative in Eq. ( |2.13| ), one finds that there exist an 
infinite number of solutions to the equation Du = 0, where m = (1 + F)u. These solutions 
are axially symmetric charge distributions, uniformly rotating in the magnetic field. But 
any finite frequency, no matter how small, excludes such solutions. 

Thus, for LU ^ Uc, Eq.( |6.4| ) turns into 
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By using asymptotical expressions for the Bessel function, one finds that for kRc S> 1 the 
interaction is screened, hke in a metal, while for kRc ^ 1 the behaviour is like in a dielectric. 
This reflects localization of electrons by a magnetic field. In reaction to a small external 
electric field, the cyclotron orbits are shifted by a small distance along the field, and also 
start drifting perpendicular to the field. Effectively, this leads to a dielectric response to an 
external electric field. 



B. Tunneling gap 

The interaction contribution to the saddle point action is J (a; — Z^) ^ J-" {uj — D) ^ J. It 
is straightforward to find that 

Sint - Sfree = — I ^ ^ J^k Jq Rc) slu^ ^ (6.6) 

Integrating over cj, one finds Sint — Sfree = where 

J m(l + F,(l-Jo^(fci?e))) ^ ^ ^ 

Since the action is linear in time, the quantity 6 is the energy of a static configuration, which 
eventually forms at t ^ l/cUc after tunneling. 

By analytic continuation of t, this contribution to the action turns K^''^^ into K^j:^^ e^"^^^ . 
Now, it is useful to integrate in Eg. ( |3.3|) , first over t, and then over dr. The integral over 
t can be easily evaluated by the residue method, since the time integration contour can be 
closed. 

It turns out that, because of the position of the poles of the contribution of 

is zero unless eV > A, while the contribution of K'^ is zero unless eV < —A. This leads to 
the result: 

IfreeiV - A) at V > A ] 

at -A < V < A ; , (6.8) 

IfreeiV + A) at V < -A 
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where I free is given by Eg. ( ^.18| ). 

Hence, the dl/dV curve has peaks at eV = ±A, and a gap of the width 2A between 
them. The apphcabihty of Eq.( |6.8| ) is restricted to the vicinity of the peaks: it correctly 
gives the width of the gap and the current near the peak, but it does not hold deep inside 
the gap, near V = 0, because our simple approach to scattering breaks down at t ^ r^. In 
this limit, the character of charge spreading changes from ballistic to diffusive, and one has 
a Coulomb anomaly problem in a diffusive system (e.g., see Ref. jS^). 



C. Qualitative discussion 

One can make a simple qualitative picture of the tunneling on the basis of Eqs.( p.l|) 
and (|6.5|) . Note that the first term in Eq. ( |6.1|) yields a solution of the non-interacting Fermi 
liquid equation, which is is localized on the cyclotron orbit in space, and in the limit uJcT ^ 1 
possesses an axial symmetry. This solution has a 5— function singularity in the momentum 
space, because at each point of the orbit the particle momentum is tangential to it. 

In the presence of interaction, the solution has an extra term, which can be obtained 
by applying to the solution for the non- interacting problem. First, since the operator 

integrates over all angles in the momentum space, the resulting contribution to 'Un(r) is 
a smooth function in the momentum space. Also, one can see that it has a spatial scale 
ao ~ ue^R^ ^ Re (at this scale the denominator in Eq.( |6.5|) becomes the order of 1). 
Therefore, the two terms of the total solution are: 

(i) a quasiparticle like 5— function singularity, localized on the classical orbit; and 

(ii) a smooth background. 



In Ref. |]T6[ the second contribution was derived hydrodynamically, while the first one was 
missed, since it corresponds to the ballistic motion of one particle and cannot be treated 
hydrodynamically. In fact, we think that the applicability of hydrodynamics for treating 
this problem is questionable, since one deals here with the collisionless limit. Nevertheless, 



our consideration justifies the picture of Ref. ITH], and corrects it only in the region r ~ i?, 
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where hydrodynamics is not apphcable. This region corresponds to the core of the "vortex" 
of Aleiner, Baranger, and Glazman. Nevertheless, it turns out that the contribution of the 
core to the gap becomes essential in the two-layer system, when the electron at one plane 
interacts with a hole on the other plane. In the one-layer system, where one should substract 
self- interaction from Eq.(|6.7D, the result of Ref. gives the correct value of the gap. 



The electrostatic energy of the core is determined by the inverse cyclotron radius (see 
Ref. [^), while the energy of the smooth background is proportional to a^^, and thus is 
much smaller. 

Let us analytically evaluate the contributions to A from 2 regions in the A;— space, by 
using the asymptotics of Jo{kRc) in the regions kRc ^ 1 and kRc ^ 1. We take the 
interaction F as the difference of the in-plane and interplane Coulomb interaction, which 
accounts for the opposite sign image in the other plane. The Fourier transform of this 
effective Coulomb interaction, with an image contribution from the other plane, is Uk = 
2Tie^{l—e~^'^)/k, and = vUk- Assuming that d <^ i?c, and Ib <^ (r^ = {2T:ve^)~^ 



is the Debye screening length, Ib = \J'hc/eB is magnetic length), one has 

A.«.«i = ^ In - (6.9) 

A,,j,»i = ^ln-^, 6.10 
'KVp Ib 

The magnetic length was chosen as an ultraviolet cutoff, since at this distance the semiclas- 
sical approximation breaks down. 

The contribution ( |6.9|) coincides with that found by a hydrodynamical approach in 
Ref. [|l^, the only difference being the effect of screening by another plane. From Eq.( |6.10D 



one can see that the region A; ^ l/r^ does give the main contribution to the gap. It can be 
interpreted as the electrostatic energy of two uniformly charged rings of charge e, radius tq, 
at distance d from each other, Ib being the effective ring thickness. These rings correspond 
to an electron and a hole uniformly spread along their cyclotron orbits. 

The contribution ( |6.1U| ) is proportional to the magnetic field. In Refs. [0,0a gap linear 



in the magnetic field was reported. Our result for the gap is in agreement with this data, 
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however, numerically it is somewhat bigger. It is natural, since the real charge distribution 
during the tunneling is more smeared than in our calculation, because Ib is not very small, 
and also because the thickness w of the 2D electron layer is comparable to the distance 
between electrons. Hence, our calculation gives only the upper bound to the gap. 

D. Effect of a parallel field 

Let us consider the effect of a magnetic field parallel to the layers. This is interesting 
because the field effect on the I — V curve reveals the spatial coherence of the tunneling. The 
parallel magnetic field does not affect the motion in the plane. However, the parallel field 
flux, "captured" between the electron and hole trajectories, will give an Aharonov-Bohm 
phase to the tunneling amplitude. Thus, the tunneling current will oscillate as the parallel 
field is swept. 

To incorporate the parallel field in the formalism, one may write its vector potential 
Az = B\\y, Ax = Ay = (i?j| is assumed to be parallel to the x— axis), and then change the 
phase of i/jr: ipR^x.y) iljji{x,y)e^^^^^'^'^^^^. (Here we restore h to show the characteristic 
scale of B\\.) Hence, the amplitude K acquires an additional factor: 

i^(r,t) = e-sfS||'^2/ir(r,t) , (6.11) 

which is just the Aharonov-Bohm phase mentioned above. For the cyclotron orbit y = 
Rc{cos9o — cos6'i), and Eg. (|3.3|) in the limit eV <^ huc yields (see Appendix B) 

I{V, S||) = J'o (^^) m B\\ = 0). (6.12) 

One may note that the tunneling current suppression by a factor oscillating as a function 
of the parallel field is similar to the magnetic field effect on the critical current in a wide 
Josephson contact. 

Note that the scale of spatial coherence Rc is much less than the scale of the charge 
spreading during the tunneling, which is given by ao ~ 1/-B^. 
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For an ideal two dimensional system, the current dependence on factors out because 
the parallel field does not affect the motion in the plane. For real wells of finite width, the 
factorization ( |6.12| ) should still be a good approximation at small B±, when the cyclotron 
radius is big compared to the well width. However, since the parallel field will squeeze the 
states in the wells, and effectively increase the barrier width, the tunneling rate may acquire 
an additional non-oscillatory suppression factor. 

VII. CONCLUSION 

The multi-dimensional bosonization technique is capable of dealing with the long wave- 
length excitations in a Fermi liquid. We applied this formalism to study tunneling into a 
Fermi liquid in a magnetic field, where one has to account for the Fermi liquid response 
to injection or removal of an electron. In the absence of the field, there is a peak of the 
tunneling conductance near zero bias. In a finite field, the peak is shifted to higher voltage, 
and a tunneling gap is formed. We calculate the gap at the field weaker than the Quantum 
Hall field, and find that it scales linearly with the field. We compare these results with the 
many-body calculation, and with experiments. 

A particular advantage of the bosonization method is that it gives an explicit relation 
with the classical dynamics of cyclotron motion, and with cyclotron orbits. We find that 
the tunneling is spatially coherent on the scale set by the cyclotron radius, and propose to 
probe this coherence by a magnetic field parallel to the barrier. The effect of the parallel 
field on the tunneling is due to an Aharonov-Bohm phase in the tunneling amplitude, which 
leads to oscillations of the current taken as a function of the parallel field. 

APPENDIX A: TUNNELING CURRENT IN LUTTINGER LIQUID 

We derive here Eq.([4.7D from Eqs.(p.3|) and ( [4.5| ). First, it is useful to integrate K{x,t) 
in Eq.( [4.5[ ) before making an analytic continuation. By use of the identity 
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dx _ r(«-i) ^^^^ 



{x^ + a^f |a|2"-i r(a) 



one gets 



a 



2(13- 



-2) (/?-2)r(/?-|) 



K{x,T)dx- 2^3/2(^1^1)2/3-3 f(^) (A2) 

(xr(x) = r(x + 1) was also used). Now, let us continue this expression from the upper 
half-axes to the upper half-plane and from the lower half-axes to the lower half-plane. For 
t > it gives: 

oo 

" {K>{x,t) - K<{x,t))dx = (A3) 

a2(/3-2) vr (/3 - 2) r (/3 - I 
- V3/2(,,)2,-3 2 - f(^^ 

Now, the integration over t can performed by using the formula 

= ^ — — , (A4) 







x" sinvrarfo; 

and the result is given by Eq.( [4.7|) 



APPENDIX B: EFFECT OF PARALLEL FIELD 



Here we derive Eg. ( 6. 12 ). We will do it for free electrons, but a generalization for a 
problem with an interaction will be straightforward, as can be seen from Eq.( |6.8|) and the 
preceding discussion. First, let us integrate over dr in Eq.( p.3|) by noticing that 



r = 2RcSm — - — , (p = — - — , (Bl) 

so that the area element rdrdip can be expressed in terms of ^ = and (f. Then one 

plugs Eqs.( |6.1I| ) and ( |5.17| ) into Eq. (|0|) . Integration over ip produces a zero-order Bessel 



function, and one arrives at 

. 2 e\to\' 7_eV 



eV [ cos —e-^^-^^' x 





TT VpUJcPF 



^ (2eB\\dRc e\ , 
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(here the interval of integration was extended to (0, oo), which accounts for the presence of 
the sum). We make a Fourier expansion of Jo(2a sin 6'/2), by using the identity 

a oo 

Jo(2asin-)= e^^'-^nW , (B3) 



and one gets 



j2 l'eBi^dRc\ 



9 U |2 oo 

/ = -^e\/V " ^ "'^ ^ , ■ (B4) 

^ hPF {eV - nu.y + ^ 

Since we are deahng with the hmit eV -C huc, only the 0— th term is essential, and one arrives 
at Eq. (|6l2D . Finally, Eq. (^J% predicts the total suppression of the tunneling current at 
certain values of B\\, given by the roots of the Bessel function Jq. However, near these points 
the current I{B\\) vanishes. This is an artifact of the approximation. If one keeps the higher 
order terms in Eg . ( P4|) , the minimal values of the current become positive. 
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